Semidilute polymer solutions at equilibrium and under shear flow 
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The properties of semidilute polymer solutions are investigated at equilibrium and under shear 
flow by mesoscale simulations, which combine molecular dynamics simulations and the multiparticle 
collision dynamics approach. In semidilute solution, intermolecular hydrodynamic and excluded 
volume interactions become increasingly important due to the presence of polymer overlap. At 
equilibrium, the dependence of the radius of gyration, the structure factor, and the zero-shear 
viscosity on the polymer concentration is determined and found to be in good agreement with 
scaling predictions. In shear flow, the polymer alignment and deformation is calculated as function 
of concentration. Shear thinning, which is related to flow alignment and finite polymer extensibility, 
is characterized by the shear viscosity and the normal stress coefBcients. 



I. INTRODUCTION 

The properties of dilute polymer solutions under shear 
flow have been studied intensively [ll-[l5j. Recent ad- 
vances in experimental single-molecule techniques even 
provide insight into the dynamics of individual polymers 
under shear flow [IHEj ■ Similarly, the dynamics of in- 
dividu al p ol y mers in a melt has been addressed exten- 
sively [16l42l| . However, we are far from a similar under- 
standing of the dynamics of semidilute polymer solutions, 
although insight into the behavior of such systems is of 
fundamental importance in a wide spectrum of systems 
ranging from biological cells, where transport appears in 
dense environments p^ . to turbulent drag reduction in 
fluid flow. Moreover, in semidilute solutions of long poly- 
mers, viscoclastic effects play an important role. Due to 
the long structural relaxation time, the internal degrees 
of freedom of a polymer cannot relax sufficiently fast 
under non-equilibrium conditions and an elastic restor- 
ing force tries to push the system towards its original 
state. Here, a deeper understanding can be achieved by 
mesoscale hydrodynamic simulations [1^ HI] . 

The dynamical behavior of dilute and semidilute poly- 
mer solutions is strongly affected or even dominated by 
hydrodynamic interactions (23l - [25| . From a theoretical 
point of view, scaling relations predicted by the Zimm 
model at infinite dilution, e.g., for the dependence of dy- 
namical quantities as viscosity and relaxation time on 
the pol yme r length, are, in general, accepted and con- 
firmed [25l427l |. However, as the concentration of the 
polymer is increased beyond the segmental overlap con- 
centration c* , where the volume occupied by polymer 
coils is equal to the total volume, the dynamics becomes 
more complex due to intermolecular excluded volume in- 
teractions [25l. [28l - l30j . For this regime, the scaling theory 
established by de Gennes describes the polymer dynam- 
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ics using the concept of 'blobs' [28|. Here a blob consists 
of g monomers and has the radius ^. A polymer chain 
comprised of Nm monomers can be regarded as composed 
of Nm/ g blobs which are hydrodynamically independent. 
Inside of a single blob, the dynamics follows the predic- 
tions of the Zimm model in dilute solution. On length 
scales larger than ^, hydrodynamic and excluded volume 
interactions are screened due to chain overlap. Thus, 
the polymer dynamics on this scale can be described by 
the Rouse model. When the concentration is further in- 
creased, the polymer dynamics is dominated by entangle- 
ment effects, which arise from physical uncrossability of 
chain segments for sufficiently long polymers. Based on 
this theory, the relaxation time and also the zero-shear 
viscosity in the semidilute regime can be scaled by using 
the concentration ratio c/c*, where c is the segment con- 
centration. Various experiments confirm the predicted 

tendencies for the relaxation time and viscosity [29l - 
However, a systematic simulation study is still lack- 
ing, even though single- chain hydrodynamic simulations 
are well established |l [13, [H-EI . The difficulty is that 
in the semidilute regime a large polymer overlap is nec- 
essary, whereas at the same time the segmental density 
has to be rather low to retain hydrodynamic interactions, 
which requires the simulation of long polymers [44| . 

The large length- and time-scale gap between the sol- 
vent and macromolecular degrees of freedom requires a 
mesoscale simulation approach in order to assess their 
structural, dynamical, and rheological properties. Here, 
we apply a hybrid simulation approach, combining molec- 
ular dynamics simulations (MD) for the polymers with 
the multiparticle collision dynamics (MFC) method de- 
scribing the solvent [H, [ll, [13, lis^tsoj . As has been 
shown, the MFC method is very well suited to study the 
non-equilibriurn properties of polymers [H, EH, |5l[ 113 , 
colloids ISjU^lljl, ^'^d other soft-matter object such as 
vesicles 55 1 and cells [5l,[53] in flow fields. 

Experiments [l|-[p, [ssl , th eoretical studies [TTI. [l5|. 
and simulations |l4l.l38l lio . [4ll . [59l - l62| of individual poly- 
mers under shear-flow conditions exhibit large deforma- 
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tions and a strong alignment of the polymers. More- 
over, a large overlap is present in a semidilute solution 
of long polymers. A typical simulation requires 10^ - 10^ 
monomers and 10^ - 10^ fluid particles. Hence, despite 
the adopted mesoscale approach, large systems can only 
be studied on a massively parallel computer architec- 
ture. Here, wc present results of large-scale simulations 
of semidilute polymer solutions under shear. The simula- 
tions were performed with our program MP^C (massively 
parallel multiparticle collision dynamics) [63| . which ex- 
hibits excellent scaling behavior on the massively parallel 
architecture of the IBM Blue Gene/P computer [gJI- For 
the MPC fluid, we find a linear increase of the speedup 
with increasing number of cores in a strong scaling bench- 
mark up to 2^^ 2", 2^6 cores for 10^, 8 x 10^6 x 10^ 
fluid particles, respectively. 

The paper is organized as follows. In Sec. II, the model 
and simulation approache are described. The equilibrium 
properties of the system are presented in Sec. HI. Sec. 
IV is devoted to the structural and conformational prop- 
erties of the system under stationary shear flow. In Sec. 
V, results are presented for the rheological properties and 
finally. Sec. VI summarizes or findings. 



into cubic cells of side length a and their relative veloc- 
ities, with respect to the center-of-mass velocity of the 
cell, are rotated around a randomly oriented axis by a 
fixed angle a, i.e., 

v,{t + h)= v,{t) + (R(a) - E)(t;,(t) - Vc^(t)), (4) 

where Vi(t) denotes the velocity of particle i at time t, 
R(a) is the rotation matrix, E is the unit matrix, and 



(5) 



is the center-of-mass velocity of the particles contained 
in the cell of particle i. N^. is the total number of solvent 
particles in that cell. 

The solvent-polymer coupling is achieved by taking the 
monomers into account in the collision step, i.e., for col- 
lision cells containing monomers, the center-of-mass ve- 
locity reads 



mNc + MNV^ 



.it) 



(6) 



II. MODEL AND PARAMETERS 

The solution consists of Np linear fiexible polymer 
chains embedded in an explicit solvent. A linear poly- 
mer is composed of Nm beads of mass M each, which are 
connected by harmonic springs. The bond potential is 
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where I is the bond length and k the spring constant. 
Inter- and intramolecular excluded-volumc interactions 
are taken into account by the repulsive, shifted, and trun- 
cated Lcnnard- Jones potential 
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where Q{x) is the Heaviside function [9(x) = for a; < 
and Q{x) — 1 for x ^ 0]. The dynamics of the chain 
monomers is determined by Newton's equations of mo- 
tion, which are integrated by the velocity Verlet algo- 
rithm with the time step hp j65j . 

The solvent is simulated by the multiparticle collision 
(MPC) dynamics method Hi [H El E!] . It is composed 
of Ns point-like particles of mass m, which interact with 
each other by a stochastic process. The algorithm con- 
sists of alternating streaming and collision steps. In the 
streaming step, the particles move ballistically and their 
positions are updated according to 



n{t + h) = n{t) + hvi(t), 



(3) 



where i = 1, . . . , and h is the time interval between 
collisions. In the collision steps, the particles are sorted 



where iV™ is the number of monomers within the con- 
sidered collision cell. To insure Galilean invariance, a 
random shift is performed at every collision step [66|. 
The collision step is a stochastic process, where mass, 
momentum and energy are conserved, which leads to the 
build up of correlations between the particles and gives 
rise to hydrodynamic interactions. 

To impose a shear fiow, for the short chains, we ap- 
ply Lees-Edwards boundary conditions [65|. A local 
Maxwellian thermostat is used to maintain the tempera- 
ture of the the fluid at the desired value 
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A parallel MPC algorithm is exploited for systems 
of long chains, which is based on a three-dimensional 
domain-decomposition approach, where particles are 
sorted onto processors according to their spatial coor- 
dinates [6^ . Here, shear flow is imposed by the opposite 
movement of two confining walls. The walls are parallel 
to the xy-plane and periodic boundary conditions are ap- 
plied in the x- and y- directions. The equations of motion 
of the solvent particles are modified by the wall interac- 
tion [g^. We impose no-slip boundary conditions by the 
bounce-back rule, i.e., the velocity of a fiuid particle is 
reverted when it hits a wall and phantom particles in a 
wall are taken into account. The same rule is applied for 
monomers when colliding with a wall [23 | . 

The simulation parameters are listed in Table ID All 
simulation are performed with the rotation angle a ~ 
130°. Length and time are scaled according to f/3 = rp/ a, 
P e {x,y,z], and t = t^ksT j {ma-^), which corresponds 
to the choice fcsT = 1, jti = 1, and a = 1, where T is 
the temperature and the Boltzmann constant. The 
collision time is /i = 0.1. The parameters yield the shear 
viscosity fj = rj/ ^raksT /a^ = 8.7. A large rotation 
angle a > 90° and a small time step h are advantages 
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to obtain high fluid viscosities, low Reynolds numbers, 
and larger Schmidt numbers. The selected values yield 
the fluid Schmidt number Sc « 14, i.e., a fluid is sim- 
ulated rather than a gas [H, HI]. Between MPC colli- 
sions, the monomer dynamics is determined by molecular 
dynamics simulations for h/hp steps, with hp ~ 0.002. 
Moreover, we choose I = a, a = a, kBT/e = 1, and 
R = no? / {ksT) = 5 X 10'^. The large spring constant k 
ensures that the mean of the bond length changes by less 
than 0.5% and the variance of the bond length distribu- 
tion by 3% only, even at the largest shear rate. 



III. EQUILIBRIUM PROPERTIES 

Before we will address polymer solutions under shear, 
the scaling behavior of equilibrium properties is dis- 
cussed, in order to determine the crossover at which a 
solution starts to follow the expected scaling laws of a 
semidilute solution. 



A. Conformational properties 



The mean square radius of gyration (Rno) in dilute 



solution obeys the scaling relation 



(7) 



with an exponent i/ w 0.59 for a good solvent [25|, |28[. 
The obtained values for (R'^q) are listed in Table U for 
various chain lengths. A fit of Eq. ([7]) to our simulation 
data obtained at the lowest concentrations (cf. Table 
1) yields the exponent t^=0.61. in good agreement with 
theory and experimental data |25l . |28| . 

As the concentration increase, the polymer coils 
start to overlap when the monomer concentration c = 
NmNp/V exceeds the value c* = Nm/Vp, with the vol- 
ume of a polymer = 4n{R'^Q)^/^ /3 and the totally 
available volume V ' 
the dependence 



25|. Scaling considerations predict 



C \ (2i'-l)/(l-3i') 



(8) 
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for the radius of gyration at concentrations c 3> c* 
[28j . This relation has been confirmed experimentally 
and by computer simulations [62l . iTOj . 

Figure [1] shows radii of gyration for various polymer 
lengths and concentrations. Our simulation results follow 
the scaling predictions. For c ^ c*, {R'^) is independent 
of polymer concentration. At c/c* « 1, the coil size starts 
to decrease and for c > c*, (Rg) ^ c""-^^^ with v — 0.61. 

As suggested by de Gennes, the coil overlap implies a 
screening of excluded volume and hydrodynamic interac- 
tions on length scales larger than the blob size ^ [2^, [2^ . 
Below this length, the swollen conformations and hydro- 
dynamic interactions are maintained. The correlation 
length ^ is independent of chain length Nm and is only 
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FIG. 1: Relative mean square radii of gyration as a function of 
the scaled concentration c/c* for the chain lengths Nm ~ 20 
(a), 40 (♦), 50 (•), and 250 (■). The solid line indicates the 
dependence of Eq. ([8} for v — 0.61. 



a function of monomer concentration at strong overlap, 
which yields the scaling relation 



-[^/(Siz-l) 



(9) 



The crossover is reflected in the polymer structure fac- 



tor 



S{q) 



N„ 



(exp[-iq(n -r^)]). 



(10) 



,1/5 



which exhibits the power-law dependence S{q) ^ 
for 1 < q{Riy/^ <. {Rl^^ll with (5=1/2 for a meh 
and 6 ~ V va. good solvent [7l| . Hence, in a semidi- 
lute solution two regimes are expected, separated by the 
correlation length ^: A good solvent behavior on length 
scales smaller than ^ and a Gaussian chain behavior on 
length scales larger than ^, i.e.. 



S{q) 



for 27r/(i?2)i/2 < ^ < 27r/^, 
for 27r/^ <q< 2tt/1, 



(11) 



In Fig. [2] polymer structure factors are shown for 
Nj-a = 250 in dilute and semidilute solutions as well as 
for ~ 50 in dilute solution. In order to obtain the 
two regimes separated by ^, the polymer chains have to 
be sufficiently long to provide not only a ratio c/c* ^ 1 
but also a low segment concentration c. As shown in the 
figure, for dilute solutions, Siq) decays with an expo- 
nent 1/^, where v = 0.61. The polymers in the semidi- 
lute regime show a crossover from the scaling behavior 
S q~^ at small q to the behavior S q~^/'^ at large 
q values. The crossover between the two regimes corre- 
sponds to q ~ 2tt/^. The values for ^ are presented in the 
inset of Fig. [2] and are found to be in excellent agreement 
with the scaling prediction © with i/ ~ 0.61. Thus, the 
scaling relation captures the concentration dependence of 
the blob size for the considered range very well. 
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TABLE I: List of simulation parameters. Lx, Ly, Lz denote the dimensions of the simulation box, 7 the shear rate, and (Nc) 
is the mean number of fluid particles in a collision cell. For Np, c, c* , and 7 the smallest and largest values used are given. 
Actual concentrations are provided in figure captions. 







Np 


Lx/ax Ly/ax L^/a 




c/r 


-3 




cjc 


7/Vfcsr/(ma2) 




20 


10 


- 200 


20x20x20 


7.05 


0.025 


- 0.5 


0.26 


0.098 - 1.96 


7.5 X 10"* 


5 


40 


10 


- 100 


20x20x20 


17.29 


0.05 - 


- 0.5 


0.13 


0.38 - 3.76 


7.5 X 10"* 


5 


50 


10 


- 512 


50x50x50 


24.51 


0.004 - 


0.205 


0.098 


0.041 - 2.08 


10"^ - 3 X 10~^ 


10 


250 


50 


- 3000 


450x75x75 


163.49 


0.0049 - 


- 0.296 


0.029 


0.17 - 10.38 


10"^ - 3 X 10"^ 


10 




FIG. 2: Structure factors of polymers of length Nm = 250 for 

the concentration ratios c/c* = 0.17 ( ), c/c* — 2.77 ( ), 

c/c* = 5.19 (- • -), and c/c* = 10.38 ( ), and of Af,„ = 50 

in the dilute regime (- •■ -). The dashed and solid straight 
lines represent power-law functions with the exponents Ijv = 
1/0.61 and 2, respectively. Inset: Dependence of the blob size 
on the concentration. 



B. Dynamics 

The polymer dynamics is dominated by hydrody- 
namic interactions in dilute solution. Theoretical re- 
sults on the concentration dependence of the relaxation 
times for small overlap concentrations are presented in 
Refs. [zl, [zll- Compared to the infinite-dilution limit, 
a term linear in the concentration is obtained, which is 
consistent with experimental data [s^. However, the ex- 
perimental data can also well be fitted by an empirical 
exponential function [soj . 

In semidilute solution, the intermolecular interactions 
between different chains become increasingly important. 
The dynamics of the polymers can be classified accord- 
ing to their intermolecular interactions as unentangled 
or entangled [29l . [30| . In the unentangled regime^ the 
monomers move according to Brownian motion in all 
three spacial directions and their dynamics can be de- 
scribed by the Rouse behavior of polymers which consist 
of blobs. Thus, the longest relaxation relaxation time 
reads as 



Tfc 



of monomers in a blob. Inside of a blob, the dynamics 
follows the Zimm behavior 



(13) 



and the longest relaxation time exhibits the concentra- 
tion dependence 



C \(2-3i^)/(3i.-l) 



(14) 



In the entangled regime^ polymers are assumed to ex- 
hibit reptation motion inside a tube caused by the pres- 
ence of neighboring chains. The monomer dynamics is 
then described by reptation theory |28| , where the longest 
relaxation time obeys the relation 



C \ (3-3f)/(3i/-l) 



(12) 



(15) 



By calculating end-to-end vector correlation functions, 
which exhibit an exponential decay, we determined the 
longest polymer relaxation times t for various concen- 
trations. The relaxation time tq at infinite dilution is 
obtained by extrapolation to zero concentration. The ob- 
tained values for tq are shown as a function of polymer 
length in the inset of Fig. |31 Their length dependence is 
well described by the power-law tq ^ , with v = 0.6, 
in accord with predictions of the Zimm model [2^ . 

Figure [3] depicts the dependence of the relaxation time 
on concentration. In the vicinity of c/c* = 1, r follows 
the scaling prediction (jl4p for an unentangled semidi- 
lute polymer solution. With increasing concentration, r 
increases faster, which we attribute to strong intermolec- 
ular interactions. Although, there are no entanglements 
in our system for c/c* < 10, the predicted dependence 
for entangled polymer melts is indicated by the solid line 
for illustration. This dependence is not reached and re- 
quires longer polymers or higher concentrations. A very 
similar dependence has been obtained experimentally in 
Ref. [soj over comparable concentration and relaxation 
time ranges. 

According to the Zimm model [MIlSjIlllzl, hydrody- 
namic interactions strongly affect the diffusive dynamics 
of polymers in solution and lead to the time dependence 



where r?, is the blob relaxation time and g the number 52 (i) = (([f"i(0~''"cm(0]~[^j(0)~'''cm(0)])^) ^ t^^'^ (16) 
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FIG. 3; Concentration dependence of longest polymer relax- 
ation times r for the polymer lengths Nm = 50 (•) and 250 
(■). The dashed line indicates the prediction Eq. p4|) and the 
continues line Eq. (|15p with u — 0.6. Inset: Polymer length 
dependence of the relaxation time at infinite dilution. The 
solid line shows the dependence tq ~ with u — 0.6. 




FIG. 4: Mean square displacements of monomers in the 
center-of-mass reference frame for the concentrations c/c* = 

0.17 (- - -), c/c* = 2.77 (- ■ -), c/c* = 5.19 ( ), c/c* = 10.38 

( ) of polymers of length Nm ~ 250. The short lines indi- 
cate the dependencies 32 ~ t^^'^ and 32 ~ t^^'^, respectively. 



of their monomer mean squared displacement in the 
center-of-mass reference frame for time scales larger than 
the Brownian time [t^ and smaller than the longest re- 
laxation time at which g2{t) saturates. In the semidilute 
regime, hydrodynamic interactions are screened for time 
scales lar ger than Tf, , the time needed to diffuse a blob 
diameter [44|. Consequently, after a time ti,, the dynam- 
ics is Rouse-like and 52 (^) ^ t^^^ Ml- Figure U displays 
(72 (i) for various concentrations for polymers of length 
= 250. 

For a dilute solution with c/c* = 0.17, 32 (0 ^ f^^^ 
in the time interval 10""^ < t/ro < 10~^ as expected. 
For t > T, 32 slowly approaches a plateau value. At 
larger concentrations, g2(t) displays a t^^^ dependence 
which turns into a t^^'^ behavior at a time Tt, which de- 
crease with increasing concentration. The concentration 
dependence of the mean squared displacement reflects the 
screening of hydrodynamic interactions in the semidilute 
regime. However, the dependence of on concentration, 
which is linked to the screening length according to 
Tb ~_£|f and $H c~'^, where 7 is predicted to be 1 
[73, [zil' 0-6 [zil, or 0.5 [Mill, respectively, cannot be 
obtained unambiguously from our simulations, because 
the different time regimes are too short. Simulations of 
longer polymers are necessary to arrive at clear and pro- 
nounced diffusion regimes. 



IV. SEMIDILUTE POLYMER SOLUTIONS IN 
SHEAR FLOW 

We now discuss the properties of polymer solutions in 
shear flow. At inflnite dilution, the flow strength is char- 
acterized by the Weissenberg number Wi = 7x0, where 
7 is the bare shear rate. For Wi ^ 1, the weak shear 
flow regime, the chains are able to undergo conforma- 



tional changes before the local strain has changed by a 
detectable amount, while for Wi ^ 1, the chains are 
driven by the flow and they are not able to relax back 
to the equilibrium conformation. This is illustrated in 
figure [5l which displays snapshots for various fiow rates. 
At small Weissenberg numbers, the polymers are only 
weakly perturbed and are close to their equilibrium con- 
formations, whereas large Wi imply large deformations 
and a strong alignment with flow. 

As pointed out in section llll Bl the polymer relaxation 
times depend on concentration. Thus, in the following, 
some properties will be characterized by the Weissenberg 
number Wic = Wir(c)/To = 7t(c). The question is, to 
what extent the influence of concentration on the poly- 
mer dynamics can be accounted for by a concentration- 
dependent Weissenberg number. As we will see, this con- 
cept applies well for all structural and dynamical prop- 
erties. 



A. Conformations 

The average shape of an individual chain in solution 
under shear is illustrated in Fig. |6]by the density distri- 
bution of monomer positions with respect to the polymer 
center of mass. At low Weissenberg numbers (Fig. (Hl^a), 
(b)), the polymers arc only weakly deformed and aligned 
with respect to the flow direction x, whereas they are 
considerably stretched and aligned in the flow direction 
and are compressed in the gradient and vorticity direc- 
tion for high shear rates (Fig. IHl^c), (d)). Figure [7] shows 
that the extent of deformation and alignment depends 
upon polymer concentration. At the same Weissenberg 
number, a larger deformation and a more pronounced 
alignment is found at higher concentrations. 
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FIG. 5: Snapshots of systems with Np = 800 polymers of 
length Nm = 250 for the Weissenberg numbers Wic = 18 
(top) and Wic = 184 (bottom). For illustration, some of the 
chains are highlighted in red. 

1. Radius of gyration 

Polymer deformation and orientation are characterized 
quantitatively by the gyration tensor 

'^^'3' = ^E^^^'.'^Ar,,^-), (17) 

2—1 

where Ar^^^ is the position of monomer i in the center- 
of-mass reference frame of the polymer. 

In Fig. [8l the relative deformation along the flow di- 
rection 

= (18) 

where G^^ = (Rg) /3 is the gyration tensor at equilibrium 
for the particular concentration, is shown for various con- 
centrations and polymer lengths. A significant polymer 
stretching appears for Wic > 1. At large shear rates, 
the stretching saturates at a maximum, which is smaller 
than the value corresponding to a fully stretched chain 
{Gxx ~ I'^^m/^'^) s-iid reflects the finite size of a polymer. 
This is consistent with experiments on DNA where 
the maximum extension is on the order of half of the con- 
tour length, and theoretical calculations [s^l . It is caused 
by the large conformational changes of polymers in shear 




x/a x/a 



FIG. 6: Monomer density distributions in the flow-gradient 
plane (a), (c) and flow-vorticity plane (b), (d) for the Weis- 
senberg numbers Wic = 1 (a), (b) and Wic = 307 (c), (d). 
The concentration is c/c* — 1.6 and the chain length is 
iVm = 50. 




x/a x/a 



FIG. 7: Monomer density distributions in the flow-gradient 
plane for the shear rate -y = 10~^ and the concentrations 
c/c* = 0.16 (a) and c/c* = 2.08 (b), which corresponds to the 
Weissenberg numbers Wic = 6.2 and Wic = 11, respectively. 
The chain length is Nm = 50. 



flow, which yields an average extension smaller than the 
contour length. Nevertheless, molecules assume totally 
stretched conformations at large Weissenberg numbers 
during their tumbling dynamics. Interestingly, a uni- 
versal dependence is obtained for SGxx as function of a 
concentration-dependent Weissenberg number Wic at a 
given polymer length, whereas in terms of Wi, polymers 
at larger concentrations exhibit a stronger stretching at 
the same Wi, as shown in the inset of Fig. |8] [12 ■ The 
latter is evident, since the longest relaxation time of a 
polymer at higher concentrations is larger and hence the 
polymer is more strongly deformed at the same shear 
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FIG. 8; Deformation ratios SGxx as function of Weissenberg 
number. Open symbols correspond to systems with Nm ~ 
50 for c/c* = 0.16 (o), c/c* = 1.6 (A), and c/c* = 2.08 
(□). Filled symbols denote results for A'^ — 250 with the 
concentrations c/c* = 0.17 (•), c/c* — 2.77 (a), c/c* = 5.19 
(♦), and c/c* — 10.38 (■). In the inset, the same data are 
shown as function of Wi. 
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FIG. 9: Relative deformations in the gradient and vorticity 
direction (inset). Open symbols correspond to systems with 
Nm. = 50 for c/c* = 2.08 (□). Filled symbols denote results 
for Nm = 250 with the concentrations c/c* — 0.17 (•), c/c* = 
1.38 (+), c/c* = 2.77 (A), c/c* = 5.19 (♦), and c/c* = 10.38 
(■)• 



rate. 

Theoretical calculations for single polymers in dilute 
solution predict the dependence SGxx = Cx^Vi^ for 
Wi < 1, where Cx is a universal constant. The renormal- 
ization group calculations of Ref. [HI yield Cx = 0.27, 
whereas a calculation based on a Gaussian phantom chain 
model yields Cx « 0.3 [u!, 113, IMIH. As shown in 
Fig. [51 the simulations confirm the quadratic dependence 
on the shear rate; SGxx is independent of chain length for 
Wi < 1 and Cx ~ 0.1. For Wi > 10, finite size effects ap- 
pear and different asymptotic values are assumed for the 
two chain lengths. We like to stress that our simulations 
are in agreement with the molecular dynamics simulation 
results of Ref. ^ and the SANS data of Refs. [H, IH. 

In the gradient and the vorticity directions, the poly- 
mers are compressed, with a smaller compression in the 
vorticity direction as shown in Fig. [5] To highlight the 
universal properties of the systems, we present the ratios 
Gpp/Cff, [P e {y,z}), where G^Oj = (i?^^) /3 is calcu- 
lated from the radius of gyration in dilute solution at 
equilibrium. At low shear rates — there is no detectable 
shear deformation by shear — polymers shrink by concen- 
tration effects for c/c* > 1 (cf. Fig. [1]). This is illus- 
trated in Fig.|n]for Wic < 10 and various concentrations. 
The dashed lines indicate the values of / (^go) from 
Fig.[TJ Evidently, the ratios Gjsp /G™ are consistent with 

the values (i?^) / (-Rgo) ^'^^ each concentration. The ra- 
tio for the shorter chain Nm = 50 and concentration 
c/c* = 2.08 is close to that of the longer chain with a sim- 
ilar concentration ratio c/c* = 2.77. This is consistent 
with the fact that (i?g) / (Rgo) is independent of chain 
length (cf. Fig. [T]). With increasing shear rate, the var- 
ious curves progressively approach a universal function, 
which decays as ~ Wi~°'^^ over the considered Wi-rangc. 
Hence, we obtain a different scaling behavior of the radius 



of gyration along the flow direction and the transverse 
directions. The reason is that a high monomer density 
is maintained along the flow direction due to polymer 
stretching, whereas the density in the transverse direc- 
tions decreases by flow-induced polymer shrinkage. 

The exponents of the power-law decay of Gyy and 
Gzz compare well with the experimental data on single 
DNA molecules Similarly, simulations (with/without 
hydrodynamic and excluded volume interactions) yield 
comparable exponents 0, |85|. However, simulations in 
Ref. for even larger Weissenberg numbers seem to 
produce an exponent closer to the theoretically expected 
value of 2 /3 [ll|,[8^ . According to theory, there is a broad 
crossover regime before the asymptotic behavior at large 
Weissenberg numbers is assumed, and the considered Wi 
fall into that crossover regime. 



2. Alignment 

The alignment of the polymers is characterized by the 
angle Xd which is the angle between the eigenvector of 
the gyration tensor with the largest eigenvalue and the 
flow direction. It is obtained from the components of the 
radius of gyration tensor via [3| 



tan(2xG) = ^ ■ (19) 

^xx ^yy 

The dependence of tan(2xG') on shear rate and concen- 
tration is shown in Fig.[TUl Again, a universal curve is ob- 
tained for the different concentrations at a given polymer 
length. Moreover, tan(2xG) seems to be independent of 
polymer length for Wic < 100, whereas we find a length 
dependence for larger Weissenberg numbers. In this high 
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FIG. 10: Dependence of tan(2xG) on shear rate. Open sym- 
bols correspond to systems with Nm ~ 50 for c/c* = 0.16 
(o) and c/c* — 2.08 (□). Filled symbols denote results 
for Nm — 250 with the concentrations c/c* = 0.17 (•), 
c/c* = 0.35 (-<), c/c* = 0.69 (T), c/c* = 2.77 (a), c/c* = 5.19 
(♦), and c/c* = 10.38 (■). The solid and dashed lines are the- 
oretical results pD. |80||. In the inset, the same data are shown 
as function of Wi. Lines are guides for the eye only. 



shear rate regime, we find tan(2xG) {Wic)~^^^- We 
like to emphasize that only the shear rate can be scaled 
in order to arrive at a universal function. The angle, or 
tan(2xG), cannot be scaled to absorb flow or polymer 
properties in an effective variable. Hence, the universal 
behavior of tlie alignment angle for various concentra- 
tions confirms that the Wcisscnbcrg number Wic is the 
correct scaling variable and that the alignment of poly- 
mers at different concentrations depends on the combi- 
nation Wic = 7T of shear rate and relaxation time only. 

The analytical description of Refs. [lli,|8Q] predicts the 
dependence 



tan(2xG) 



\LWi* 



1/3 



(20) 



for semiflexible polymers in dilute solution in the limit 
Wi* oo. Here, we introduce the Weissenberg number 
Wi* = 7rth for the theoretical result, because the relax- 
ation times from theory and simulation might not be the 
same; L is the length and Ip the persistence length of 
the polymer. The analytical result describes the simula- 
tion data well at large shear rates, when the Weissenberg 
number of the theoretical model is set to Wi* = Wic/2. 
To compare the predicted length dependence with that 
of the simulation, we apply the relation (i?go) = 2lpL 
to obtain a persistence length, with {R^q) the polymer 
mean square end-to-end distance in dilute solution at 
equilibrium, which yields Ip/L w 0.025 for Nm = 50 and 
Ip/L w 0.008 for Nm = 250. With these values, the ra- 
tio of tan(2x(3) of the polymer of length L = 50a and 
L = 250a is 1.5. This compares well with the factor 1.33 
following from the simulation results, which suggests that 
excluded volume interactions are of minor importance for 
intermediate flow rates. 



FIG. 11: Comparison of experimental and simulation data 
for the average orientation angle xg- The concentrations are 
c = 0.016/r^ (•), c = 0.205/r^ (■) for our simulations (with 
iVm = 50), and O.llSff// (o), 1.094^// (□) for the experiments 



In the limit Wic — > 0, theory predicts tan(2xG) ^ 
Wi~^. The simulation data do not show this dependence 
on the considered range of Weissenberg numbers, which 
might be due to excluded volume interactions. 

The inset of Fig. [TU] displays a strong dependence of 
Xg on concentration. The higher the concentration, the 
more the chains are orientated along the flow direction. 
Such a concentration effect has also been reported in 
light scattering experiments [s^, where dilute polymer 
solutions are considered. A comparison of the experi- 
mental data with the simulation results is presented in 
Fig. [TTJ Two data sets are presented, a dilute solution, 
with the concentration 0.113g/l, and a semidilute solu- 
tion with the approximately ten times higher concentra- 
tion l.OMg/l, both taken from Fig. 8 of Ref. [H]. Evi- 
dently, the experimental data fit well with our simulation 
results. Moreover, both, experiments and simulations, 
yield a shift of the curves for the higher concentrations 
to smaller Weissenberg numbers. In Ref. [s^, a factor 
Pe ~ ['7]7, where [77] is the intrinsic viscosity, is used to 
present the data. This quantity is proportional to Wi, 
since [rj] is proportional to the longest relaxation time r; 
however, f3e and Wi arc not identical. No adjustment pa- 
rameter is used in Fig. [Til which suggests that the ratio 
/3e/Wi is close to unity. 



V. RHEOLOGY 

A. Shear viscosity 

Under shear flow, the viscosity 77(7) is obtained from 
the relation 



'7(7) 



(21) 



9 



where Uxy is the shear stress 
(J, 



|88[. In our simulations, 
xy is calculated using the virial formulation of the stress 



tensor [65|,|68!,|89|. For sufficiently weak flow, the polymer 
solution is in the Newtonian regime, i.e., Gxy ~ 7 and the 
viscosity is independent of shear rate. Thus, the viscosity 
obtained in this low shear rate regime is equal to the zero- 
shear viscosity denoted by 770 • The latter depends on the 
monomer concentration, which is often presented in the 
form 



(22) 



where ?2s is the solvent viscosity and the Huggins con- 
stant [3ll[83. 

We determine the intrinsic viscosity by a linear extrap- 
olation to zero concentration of both, (770 ~ Vs)/cris and 
the inherent specific viscosity [ln(ryo/?ys)]/c. The common 
intercept of these two functions gives [rj] [3l|, [oO] ■ 

The intrinsic viscosity is proportional to Rg/Nm [111 
and is therefore proportional to the inverse of the over- 
lap concentration c* |9l|. We find ['n]c* « 0.9 and 
[77] c* « 1 for the polymer of length N,n = 40 and 
Nm = 50, respectively, which means that the propor- 
tionality coefficient is close to unity for the considered 
model systems. 

The Einstein relation 



77 = r;, (1 + 2.50) , 



(23) 



where (j) is the volume fraction, captures the concentra- 
tion dependence of hard sphere suspensions for (f> <^ 1. 
This relation should also apply for dilute polymer solu- 
tions, when the hydrodynamic radius Rh is used to define 
the volume fraction, i.e., (j) = {4:Tr/3)R^Np/V . Equa- 
tions ^ and (HSl) are consistent if [r]]c* = 2.5{RH/Rgf- 
Since [i]]c* « 1, as explained above, consistency requires 
Rg/Rn ~ 1-36. From our simulations, we find the hydro- 
dynamic radii Rh ~ 3.8Z and Rh ~ 9.4Z for the polymers 
of length N,n = 50 and N,n = 250, respectively, which 
yields the ratios 1.3 and 1.36. These values are very close 
to the value necessary to match the Einstein relation. 
The ratios are somewhat smaller than the asymptotic 
value Rg/Rn « 1.59 for N„i 00 obtained in Ref. [82], 
which is a consequence of the fact that we consider in- 
sufficiently long chains. 

The term fc// ([ryjc)^ depends on hydrodynamic interac- 
tions. The value of the Huggins constant kn of flexible 
polymers is in the range of 0.2 — 0.8 and depends on sol- 
vent quality [3l|,[9^. In good solvent, typically the value 
0.3 is found experimentally [90j . Expressing Eq. (|22l) in 
terms of the dimensionless parameter [77] c as 



VR 



1 + fcffWc- 



(24) 



which is denoted as relative viscosity, allows us to deter- 
mine kn ■ The inset of Fig. [121 shows 777^ — 1 as function of 
[77] c for polymers of length iV,„ = 40 and = 50. The 
slope of the solid line is kjj = 0.35, in close agreement 
with experiments [3l|, |9^ . 




FIG. 12: Dependence of the zero shear viscosity on the scaled 
concentration c/c* for the polymer lengths Nm = 40 (a), 
iVm = 50 (•), and Nm = 250 (■). The solid line indicates the 
power-law (c/c*)^/''^'^^^' with u = 0.6. In the inset, rjn — 1, 
Eq (|24p . is shown as function of [ri]c for N,n — 40 (A) and 
Nm = 50 (•); the slope of the solid line is 0.35, which corre- 
sponds to the Huggins constant of polymers in good solvent. 
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FIG. 13: Dependence of the polymer contribution to shear 
viscosity on shear rate. Open symbols correspond to systems 
with Nm = 50 for c/c* = 0.16 (o), c/c* = 1.6 (O), and c/c* = 
2.08 (□). Filled symbols denote results for Nm = 250 with the 
concentration c/c* = 0.17 (•), c/c* — 0.35 {■<), c/c* = 1.38 
(►), c/c* = 2.77 (A), c/c* = 5.19 (♦), and c/c* = 10.38 (■). 



For semidilute unentangled polymer solutions, the vis- 
cosity is proportional to the number of blobs per chain 
and can be expressed by the scaling relation [2^, [2^ |3l|, 

lai 



1/(31.-1) 



(25) 



Figure [T^ displays zero-shear viscosities as function of 
concentration for various polymer lengths. For c/c* > 3, 
the data are close to the power-law of Eq. pS]) . 

At sufficiently large shear rates, the polymers are 
aligned and deformed, which implies shear thinning [TTI . 
[T^ . Isol . [ssl . Figure [T3] shows the polymer contribution 



10 



■qP to the shear viscosity. Similar to the ahgnment angle, 
the viscosity is a universal function of the Weissenberg 
number Wic and shows a weak dependence on polymer 
length. It is independent of shear rate for Wic ^ 1, 
decrease approximately as WiJ°'^ for 1 < Wic < 10^, 
and Wi~"'^^ for higher shear rates. This behavior is con- 
sistent with other simulation results ^ |H, H^, [s^ [92] ■ 
However, an even stronger decay of the viscosity is ob- 
served in simulations at larger shear rates in Refs. (9l.[l3|. 
Experiments of dilute polymer solutions reported expo- 
nents ranging from —0.4 to —0.85 [2, [s^. Theoretical 
calculations for dumbbells and finite extendable poly- 
mers predict the dependence T]p Wi~^^'^ in the limit 
Wi ^ oo [nl [H dfl m, m. The differences in the ob- 
served behavior can be explained by a broad crossover 
regime before the asymptotic behavior is reached. 

The ratio of the viscosities of the two lengths is approx- 
imately 1.33 for the large Weissenberg- number regime, 
as for the alignment angle, which compares well with the 
theoretically predicted length dependence in Eq. (^0)) (cf. 
Sec. IV A. 2). 

B. Normal stress coefficient 

The concentration and shear rate dependencies of the 
first and second normal stress difference [H, |9^ 

^-l =(tT,, - <Jyy)/j^, (26) 
*2 H'^yy - <Jzz)/i^ (27) 

are displayed in Fig. [T31 Within the accuracy of the 
simulations, the ratio where VP]* is the stress dif- 

ference at zero shear rate, is an universal function of Wic 
for various concentrations and decreases as ~ j-^/^ 
for large shear rates. This is consistent with analyti- 
cal calculations ISOl. ISSl 1931 . various computer simula- 
tions d, m, lU |94| . and experiments [1, Q for 
dilute solutions. Similar to the viscosity, the decrease is 
related to the finite polymer extensibility. Both, hydro- 
dynamic and excluded volume interactions contribute to 
^1, which is shown in Ref. [1^ for single polymers. Here, 
wc find the power law 

*?-(^)'^ (28) 

as shown in the inset of Fig. fMla). Hence, the first nor- 
mal stress difference exhibits a significant dependence on 
excluded volume interactions. At large concentrations, 
might saturate; at least, we cannot exclude such a 
saturation at the upper limit of the considered concen- 
tration range. 

Second normal stress differences are presented in 
Fig. [HT b) for various concentrations. At low concen- 
trations, their values are much smaller than those of 
and hence cannot be calculated within the same accu- 
racy, and the values are difficult to find. Wc therefore 
present the simulations results for 4'2 directly rather than 




FIG. 14: First and second normal stress coefficients (a) 
and \I'2 (b) for polymers of length Nm = 50 and the concen- 
trations c/c* = 0.16 (o), c/c* = 0.41 (x), c/c* = 0.81 (+), 
c/c* = 1.63 (O) and c/c* = 2.08 (□). The solid lines indicate 
the power-law decay 'ifi ~ Wi^T*^'*. Inset in (a): Concen- 
tration dependence of the zero-shear-rate first normal stress 
coefficient Inset in (b): Ratio of \E'2/*I'i for the various 
concentrations. 



in scaled form. Similar to ^'i, the second normal stress 
difference decreases as 4*2 ~ with increasing shear 

rate. Again, excluded volume and hydrodynamic inter- 
actions contribute to ^'2 [H, [H, (H] . The ratio '^2/'^! 
is concentration dependent, as shown in Fig. [HT b). At 
small Wic and large concentrations, the ratio is close to 
unity, decreases with increasing shear rate and assumes 
a constant value above a certain Wic, which seems to de- 
pend on concentration. The plateau value itself increases 
with increasing concentration. A similar plateau has 
been found in simulations of dilute solutions in Ref. [1^ . 
The concentration dependence of the plateau value sug- 
gests that excluded volume interactions determine the 
behavior of the normal stress differences. If hydrody- 
namic interactions would be dominant, we would expect 
a decrease of the plateau with increasing concentration 
due to screening of hydrodynamic interactions by poly- 
mer overlap. 
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VI. CONCLUSIONS 

We have calculated conformational, dynamical, and 
rheological properties of polymers in dilute and scmidi- 
lute solution under shear flow by mesoscale hydrody- 
namic simulations. At equilibrium, our simulations con- 
firm the scaling predictions for the concentration depen- 
dence of the radius of gyration and the longest relaxation 
time. Moreover, we find signatures for the screening of 
excluded volume interactions in the static structure fac- 
tor. 

In shear flow, the polymers exhibit deformation — the 
polymer is stretched in flow direction and shrinks in the 
transverse directions — and alignment, which depend on 
shear rate and concentration. As one of the the main 
results of the paper, we have shown that the relative 
deformation SGxx in the flow direction, the alignment 
tan(2xG)j and the viscosity rj/rjQ are universal func- 
tions of the concentration-dependent Weissenberg num- 
ber Wic = 7^(0) [93, • This is surprising because r 
increases rapidly with increasing concentration. More- 
over, it indicates that the dynamics under shear flow is 
still governed by the relaxation time at equilibrium de- 
spite the anisotropic deformation of a polymer. This is 
not evident a priori, as expressed by the scaling behav- 
ior of the radius of gyration tensor components Gyy and 
Gzz- Here, we find a concentration-independent scaling 
behavior at large Wic only when these values are scaled 



by their equilibrium values in dilute solution. Hence, the 
deformations transverse to the flow directions seem to 
exhibit a scaling behavior corresponding to a dilute so- 
lution, however, with the relaxation time of the concen- 
trated system. 

In addition, the zero-shear viscosity obeys the scal- 
ing predictions with respect to the concentration depen- 
dence. Moreover, for the first time, we show by simula- 
tions that the Huggins constant is equal to kn ~ 0.35 
for a flexible polymer in good solvent, which is in close 
agreement with experimental results [90| . 

Finally, we find a strong concentration dependence of 
the normal stress differences. Their ratio shows that 
intermolecular excluded volume interactions determine 
their behavior at all shear rates. 

Our simulations reveal a complex interplay between 
shear rate, deformation, and intramolecular excluded vol- 
ume interactions, which is difhcult to grasp by analytical 
theory. 
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